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Abstract 

In ref. |l[ it is claimed to have proved that Lagrangian theories with a nonlocality of finite 
extent are necessarily unstable. In this short note we show that this conclusion is false. 

1 Introduction 

In ref.[P a canonical formalism for nonlocal Lagrangians with a nonlocality of finite extent is 
established. It is compared with Ostrogradski formalism |Q for local Lagrangians which depend 
on a finite number of derivatives of coordinates. One of its central conclusions is that Lagrangian 
systems with a nonlocality of finite extent have no ". . . possible phenomenological role [. . . ] . They 
have inherited the full Ostrogradskian instability . . . " . 
The aim of the present note is: 

(i) to point out some defects in ref. jl| (section ||) concerning the application of the variational 
principle that underlies the derivation of the nonlocal equation (2) from a Lagrangian; 

(ii) to stress the importance of the functional space where the variational problem is developed, 
this is also the functional space where the solutions must be searched (section ||) and 

(iii) to illustrate by two simple counterexamples (section ||) that: (A) Lagrangian systems con- 
taining derivatives of a higher order than first are not necessarily unstable and (B) nonlocality 
of finite extent does not inevitably lead to instability. 

2 The nonlocal action principle 

Although the canonical formalism set up in ref. 0| is derived on a general ground, it is basically 
illustrated by the simple nonlocal Lagrangian system 

L[q]{t) - + y) - ^^^'^(O + A) (1) 
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and the equation of motion for this Lagrangian is written as: 
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-m \ q{t) + -iu^qit + A) + -iu''q{t - A) f> = 
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It must be noticed that the latter equation as it reads does not properly correspond to a 
standard action principle. Indeed, the action integral whose variation would be the left hand side 
of (0) is: 
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where t is the same both in the numerator and the denominator. However, the Euler-Lagrange 
equation that follows from the action principle SS = is 



Sq{t') 



0, 



which is much more restrictive than (j^. 
Moreover, an equation like 



m \ q{t) + ^LO^q{t + A) + ^Lo\{t - A) \ = , 
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valid for — oo < i < oo, cannot be derived from an action integral like (^), extending over a finite 
interval. Indeed, the variation of the action (0) is: 
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The extremal condition 5S = Q then leads to the boundary conditions 5q{t + ■§■) = 5<l{t 
and to the equations of motion: 
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(d) t + ^ <T <t + A 
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which has the only solution g(r) = , for i — A < r < i + A, as it ineluctably follows from 
sequentially exploiting (d), (b), (a) and (c). 

Furthermore, if we alternatively try with an action extended over a larger interval: 

S = [ dTL[q]iT) 
Jo 

the Eulcr-Lagrange equations are: 
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[g(T + A)+g(r-A)] -0 
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where the domains {i) to {v) respectively correspond to: < r < 
r<T<T+fandT+f<r<r+A. Equation (|) only looks Ukc (|0| 
The conditions (||.i) and (|[v) then yield: 
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in the interval A < r < T. 
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that act as constraints on the possible solutions of (^ii), (^-iii) and (||-iv). As a consequence, 
equations (|^) can be reduced to an ordinary differential equation, whose order depends on the 
number of times that the elementary length A fits into [0,T]. 

We have thus illustrated the important role played by the integration bounds in the nonlocal 
action (|^) as far as the Euler-Lagrange equations are concerned. The integration bounds in the 
action and the problems associated to them are commonly overlooked in theoretical physics liter- 
ature because, in standard local cases no trouble is usually entailed by proceeding in this manner. 
Nonlocal cases are however a new ground where nothing can be taken for granted. 

For a local action, the bounds of the integral also determine the functional Banach space where 
the variational calculus is meaningful 0, e. g., the space C^([a, b]) for an action integral extending 
over [a, 6]. This is also the space where the solutions to the Euler-Lagrange equations have to be 
sought. 

A way to derive equation (H), for t extending from —oo to cx), from an action principle could 
consist in taking the integral over the whole IR: 



5 = 
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but then two additional difficulties arise: on the one hand, the action S does not converge anymore 
for all q G C^(R) and, on the other, C^(IR) is not a Banach space. (The variational calculus should 
be then approached in terms of Frechet spaces To my knowledge, it remains an open 

problem to establish the appropriate mathematical framework where a nonlocal equation like (H) 
can be derived from an action integral like (^. This results in a lack of preciseness in the definition 
of the functional space where the nonlocal equation has to be solved. 
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3 The stability problem 

Leaving aside the difficulties just mentioned, suppose that, for some physical reasons whatsoever, 
we are only interested on the solutions of in the Banach space 

B= {qeC^{R); \q{t)\, \q{t)\ and \q{t)\ are bounded} . 

The general solution of (||) is thus 

q{t) = {Aie''"' + A^e-**) (10) 

where ±fc; are the real solutions [| of 

h{k) = - Lu^ cos{kA) = (11) 

and A* is the complex conjugate of Ai, to ensure that q{t) e R. 

Notice that the number of real roots of is finite. A look at figure |l| is enough to get 
convinced that they can be indexed so that 

kj < ki if j<i and l^l,2,...N. 

N can be either odd and then all roots are simple, or even, in which case ±ki\r are both double. 
It should also be remarked that the greater is A, the denser is the wiggling in the graphics (figure 
|l|). Therefore, N increases with A. 

h(k) = k^- (D cosf/tAj 




Figure 1: h{k) = has a finite number of real roots, and the sign of the derivative h'{ki) at each 
root is alternating. 

The space of solutions of (^) in B can be hence coordinated by 2N < oo real parameters, 
namely, the real and imaginary parts of Ai, that can be put in correspondence with the initial 

, , ' . (2Af-l) 

data: go, 90, ■ ■ 

The solutions of (^) in B are stable because a small change in the initial data SqQ^^ , results in 
a small change in the complex parameters: 6Ai . Indeed, from the linearity of eq. (^) and from 

complex value of fc; would result in an exponential growth cither at +co or — oo and then q ^ B 
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the general solution ( p^ it follows that: 



N 

J2[5Ai {ikir + 5A*{-ikir] = 5qi''\ 
1=1 

a = 0, 1, ... , 2N — 1 , which can be inverted to obtain SAi as a linear function of Sqq^K Therefore, 
there exists K > such that \SAi\ < K\\Sqo\\ , where \\Sqo\\ = sup{|(5g^"'' | ; a = 0, 1, . . . 27V - 1} . 
The deviation from q(t) evolves with time as: 



\6q{t)\ = 



N 



1=1 



N 



<J2^\SA\ <2NK\\Sqo\ 



1=1 



which proves the stability f\ of the solutions of (||) in the space B. 

Notwithstanding, if we now have a look at the Hamiltonian (equation (48) in ref. Q 

H{t) = ^mq\t) + ^mu;\{t)q{t + A) - 

—muP' I dsq{t + s)q{t + s — A) , 



2 

on substituting the general solution (p^), we obtain that 



N 



H{t) =2mJ2ki9{ki)AiA* (12) 
1=1 

with ^ 

g{k) = fc + Y A sin(fcA) (13) 

As expected, H{t) is an integral of motion, but it has not a definite sign. Indeed, notice that 
g{k) = h'{k)/2 alternates sign at each root fc; [see eq. ( pi] ) and figure Therefore, g{ki) is 
positive or negative depending on whether I is even or odd, respectively (moreover, g{ki) = Q ii ki 
is a double root). 



4 Two simple counterexamples 

4.1 The so-called Ostrogradskian instability 

In ref. it is proved that the Hamiltonian formalism for a nonlocal Lagrangian can be obtained 
as a limit case for — > oo of the Ostrogradski formalism ||] for a Lagrangian that depends on the 
derivatives of the coordinates up to order N . ^ For iV > 1, the Ostrogradski Hamiltonian is linear 
on all the canonical momenta but one, namely. Pi, ... P/v_i, therefore it has not a definite sign. 

^In the sense of Liapounov, see |y 

similar result was also obtained in Jaen, X., Jauregui, R., Llosa, J. and Molina, A., Phys. Rev. D 36, 2385 

(1987) 
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The fact that the energy is not bounded from below is then argued to conclude that the solutions 
of the equations of motion are ineludibly unstable. This is what is called the Ostrogradskian 
instability. It is also shown in Q that this drawback also holds in the limit N oo. 

Actually what has been proved there is only that the energy cannot be taken as a Liapunov 
function |^ to conclude the stability of the equations of motion derived from an A^*'' order La- 
grangian {N > 1). However, the fact that a sufficient condition of stability is not met does not 
imply instability. Let us consider the following simple counterexample: 



Liq,<i,<i) = -^f 



2 ^ 



(14) 



where B and C are two parameters which we shall later tunc in order to get stability. 

According to Ostrogradski theory, the canonical coordinates and momenta are [in the notation 
of ref. [0, eqs. (6-7)] 

Qi=q, Q2 = q Fi=Sg-g(™), P2 = q 
and the Hamiltonian is (eq.(9) in (l)): 

H = -P^ 



PiQ2-\bqI-\cq\ 



(15) 



Introducing 



X = 



Qi ^ 

Q2 
Pi 



and 



IM 



/ 



c 



\ P2 J \ B -I J 

the Hamilton equations for (|l5|) can be then written as the linear system: 

d 



dt 



X = IMA , 



the stability of whose solutions depends on the real part of the roots of the characteristic polynomial 
Pm(A) = det(IM - AII4), that is: 



Ib±VB^-4C 



If the parameters are tuned so that B < and < C < then all roots are imaginary and 

the system is stable Q. The latter is not an obstacle for the fact that the Hamiltonian has not a 
definite sign. 



4.2 A case of finite extent nonlocality 

In the next example, the boundaries of the action integral are finite. Consider the nonlocal action 

S[q]^J^dtL[q]{t), with 



L[qm^y^t)-lc.'q^it) 



q{t) / dt'Git,t')q{t') 



(16) 
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where, for {t,t') G [0,T]2, 

SYauj[T - t) sm.ujt' e{t ~ t')] (17) 



Git.t') = r-^—[sinLo(T-t')s\nujte(t'-t) + 



and is the sohition of: 

dfG{t, t') + uj^G{t, t') = 6{t - t') (18) 

for the boundary conditions: G(0, t') = G{t, T) = 0. 

The variation 5S = i) with the boundary conditions 5q{Q) — 5q{T) — leads to the equations 
of motion ^ 

q{t)+LO^q{t)-LO^ ( dt'G{t,t')q{t') ^0 (19) 

The solutions q{t) must be sought in the Banach space C^([0,T]). 

Differentiating twice (^9|) and taking (|l^ into account, we arrive at: 

+ 2cj2 ^- ^ (20) 

Hence, the solutions of ( |l9| ) must be among the general solution of (|20|): 

q{t) = Ae^"* + A*e-'°'' + Dt + E (21) 

with a = wV2. The parameters A, A*, D and £^ must fulfill the following constraints 

£1 = and E = A + A* 

which result from substituting (^) into (|l9|). 
The general solution of (|9|) is therefore: 

q{t) = A (e^"* + 1) + ^* (e-'"* + 1) (22) 

The phase space for our system is thus two-dimensional, and every solution is determined by the 
initial values qo and qo: 

qo = 2{A + A*) and qo = ia{A - A*) 

By direct inspection of (|2^), we see that the solutions of equation are stable, although the 
latter is derived from a Lagrangian with a nonlocality of finite extent. That is, for any e > there 
exists p > such that \6qo\ + \6qo\ < p implies that |^g(i)l + 1^9(01 < which proves the 

stability. 



5 Conclusion 

We have intended to stress the crucial importance of clearly precising the Banach space where the 
variational principle for a nonlocal Lagrangian is formulated. This degree of precision is usually 
obviated in theoretical physics (i. e., for local Lagrangians) without any major problem. However, 
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such non rigourous way of proceeding cannot be extrapolated to systems with a new complexity. 
The relevance of the above mentioned Banach space is twofold: (i) it is where the solutions of the 
equations of motion must be sought and (ii) it is the function space where path-integrals are to be 
calculated in an eventual quantization of the system. 

We have also analised the stability of the equations of motion for a Lagrangian system presenting 
a nonlocality of finite extent. We have shown that the choice of the Banach space where the 
variational principle is meaningfully formulated is crucial to decide the stability or unstability of 
the system. Furthermore, we have seen that a system can be stable in spite of the fact that the 
Hamiltonian does not have a minimum. 

Finally, we have shown by a counterexample that higher order Lagrangian systems are not 
necessarily unstable. The fact that a sufficient condition for stability is not fulfilled does not imply 
instability. 
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